Optimal Orthogonal Graph Drawing with Convex Bend Costs* 



Thomas Blasius, Ignaz Rutter, Dorothea Wagner 

Karlsruhe Institute of Technology (KIT) 
f irstname . lastnameOkit . edu 



Abstract 

Traditionally, the quality of orthogonal planar drawings is quantified by cither the total 
number of bends, or the maximum number of bends per edge. However, this neglects that in 
typical applications, edges have varying importance. Moreover, as bend minimization over all 
planar embeddings is TV'P-hard, most approaches focus on a fixed planar embedding. 

We consider the problem OptimalFlexDraw that is defined as follows. Given a planar 
graph G on n vertices with maximum degree 4 and for each edge e a cost function cost e : No — > K 
defining costs depending on the number of bends on e, compute an orthogonal drawing of G of 
minimum cost. Note that this optimizes over all planar embeddings of the input graphs, and 
the cost functions allow fine-grained control on the bends of edges. 

In this generality OptimalFlexDraw is AfV-hard. We show that it can be solved efficiently 
if 1) the cost function of each edge is convex and 2) the first bend on each edge does not cause any 
cost (which is a condition similar to the positive flexibility for the decision problem FlexDraw). 
Moreover, we show the existence of an optimal solution with at most three bends per edge 
except for a single edge per block (maximal biconnected component) with up to four bends. For 
biconnected graphs we obtain a running time of 0(n ■ Tfl ow (n)), where Tg ow (n) denotes the time 
necessary to compute a minimum-cost flow in a planar flow network with multiple sources and 
sinks. For connected graphs that are not biconnected we need an additional factor of 0(n). 

1 Introduction 

Orthogonal graph drawing is one of the most important techniques for the human-readable visu- 
alization of complex data. Its aesthetic appeal derives from its simplicity and straightforwardness. 
Since edges are required to be straight orthogonal lines — which automatically yields good angular 
resolution and short links — the human eye may easily adapt to the flow of an edge. The readability 
of orthogonal drawings can be further enhanced in the absence of crossings, that is if the underlying 
data exhibits planar structure. Unfortunately, not all planar graphs have an orthogonal drawing in 
which each edge may be represented by a straight horizontal or vertical line. In order to be able to 
visualize all planar graphs nonetheless, we allow edges to have bends. Since bends obfuscate the 
readability of orthogonal drawings, however, we are interested in minimizing the number of bends 
on the edges. 

In this paper we consider the problem OptimalFlexDraw whose input consists of a planar 
graph G with maximum degree 4 and for each edge e a cost function cost e : No — > M defining costs 
depending on the number of bends on e. We seek an orthogonal drawing of G with minimum cost. 
Garg and Tamassia [2] show that it is A^P-hard to decide whether a 4-planar graph admits an 
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orthogonal drawing without any bends. Note that this directly implies that OptimalFlexDraw 
is A^P-hard in general. For a special case, namely planar graphs with maximum degree 3 and 
series-parallel graphs, Di Battista et al. [3j give an algorithm minimizing the total number of bends 
optimizing over all planar embeddings. They introduce the concept of spirality that is similar to 
the rotation we use (see Section 2.3 for a definition). Blasius et al. [2] show that the existence 
of a planar 1-bend drawing can be tested efficiently. More generally, they consider the problem 
FlexDraw, where each edge has a flexibility specifying its allowed number of bends. For the case 
that all flexibilities are positive, they give a polynomial-time algorithm for testing the existence of 
a valid drawing. 

As minimizing the number of bends for 4-planar orthogonal drawings is A^P-hard, many results 
use the topology-shape-metrics approach initially fixing the planar embedding. Tamassia |15| de- 
scribes a flow network for minimizing the number of bends. This flow network can be easily adapted 
to also solve OptimalFlexDraw even for the case where the first bend may cause cost, however, 
the planar embedding has to be fixed in advanced. Biedl and Kant |T] show that every plane graph 
can be embedded with at most two bends per edge except for the octahedron. Morgana et al. |12j 
give a characterization of plane graphs that have an orthogonal drawing with at most one bend per 
edge. Tayu et al. [17] show that every series-parallel graph can be drawn with at most one bend 
per edge. All these results and the algorithm we present here have the requirement of maximum 
degree 4 in common. Although this is a strong restriction it is important to consider this case since 
algorithms dealing with higher-degree vertices (drawing them as boxes instead of single points) rely 
on algorithms for graphs with maximum degree 4 |16| El [TT] . 

Even though fixing an embedding allows to efficiently minimize the total number of bends (with 
this embedding), this neglects that the choice of a planar embedding may have a huge impact on 
the number of bends in the resulting drawing. The result by Blasius et al. [2] concerning the 
problem FlexDraw takes this into account and additionally allows the user to control the final 
drawing, for example by allowing few bends on important edges. However, if such a drawing does 
not exist, the algorithm solving FlexDraw does not create a drawing at all and thus it cannot 
be used in a practical application. Thus, the problem OptimalFlexDraw, which generalizes the 
corresponding optimization problem, is of higher practical interest, as it allows the user to take 
control of the properties of the final drawing within the set of feasible drawings. Moreover, it allows 
a more fine-grained control of the resulting drawing by assigning high costs to bends on important 
edges. 



Contribution and Outline. Our main result is the first polynomial-time bend-optimization 
algorithm for general 4-planar graphs optimizing over all embeddings. Previous work considers 
only restricted graph classes and unit costs. We solve OptimalFlexDraw if 1) all cost functions 
are convex and 2) the first bend is for free. We note that convexity is indeed quite natural, and 
that without condition 2) OptimalFlexDraw is ATP-hard, as it could be used to minimize the 
total number of bends over all embeddings, which is known to be AAP-hard |9j. 

In particular, our algorithm allows to efficiently minimize the total number of bends over all 
planar embeddings, where one bend per edge is free. Note that this is an optimization version 
of FlexDraw where each edges has flexibility 1, as a drawing with cost exists if and only 
if FlexDraw has a valid solution. Moreover, as it is known that every 4-planar graph has an 
orthogonal representation with at most two bends per edge [T] , our result can also be used to create 
such a drawing minimizing the number of edges having two bends by setting the costs for three or 
more bends to oo. 

To derive the algorithm for OptimalFlexDraw, we show the existence of an optimal solution 
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Figure 1: (a) Two parallel edges, the thin has one bend for free, every additional bend costs 1, the 
thick edge has two bends for free, every additional bend costs 2. Whether embedding E\ or £2 is 
better depends on the number of bends. The minimum (marked by gray boxes) yields a non-convex 
cost function, (b) The non-convexity in (a) does not rely on multiple edges, the thick edge could 
be replaced by the shown gadget where each edge of the gadget has one bend for free and every 
additional bend costs 2. (c) This example has a non-convex cost function even if every edge has 
one bend for free and each additional bend costs 1. 



with at most three bends per edge except for a single edge per block with up to four bends, 
confirming a conjecture of Rutter |14| . 

Our strategy for solving OptimalFlexDraw for biconnected graphs optimizing over all planar 
embedding is the following. We use dynamic programming on the SPQR-tree of the graph, which 
is a data structure representing all planar embeddings of a biconnected graph. Every node in 
the SPQR-tree corresponds to a split component and we compute cost functions for these split 
components determining the cost depending on how strongly the split component is bent. We 
compute such a cost function from the cost functions of the children using a flow network similar 
to the one described by Tamassia [TS]. As computing flows with minimum cost is A^P-hard for 
non-convex costs we need to ensure that not only the cost functions of the edges but also the 
cost functions of the split components we compute are convex. However, this is not true at all, 
see Figure [T] for an example. This is not even true if every edge can have a single bend for free 
and then has to pay cost 1 for every additional bend, see Figure [ljc). To solve this problem, we 
essentially show that it is sufficient to compute the cost functions on the small interval [0,3]. We 
can then show that the cost functions we compute are always convex on this interval. 

We start with some preliminaries in Section[2] Afterwards, we first consider the decision problem 
FlexDraw for the case that the planar embedding is fixed in Section [3] In this restricted setting 
we are able to prove the existence of valid drawings with special properties. Blasius et al. [2] show 
that "rigid" graphs do not exist in this setting in the sense that a drawing that is bent strongly can 
be unwound under the assumption that the flexibility of every edge is at least 1 . In other words this 
shows that graphs with positive flexibility behave similar to single edges with positive flexibility. 
We present a more elegant proof yielding a stronger result that can then be used to reduce the 
number of bends of every edge down to three (at least for biconnected graphs and except for a single 
edge on the outer face). In Section [4] we extend the term "bends", originally defined for edges, to 
split components and show that in a biconnected graph the split components corresponding to the 
nodes in its SPQR-tree can be assumed to have only up to three bends. In Section [5] we show that 
these results for the decision problem FlexDraw can be extended to the optimization problem 
OptimalFlexDraw. With this result we are able to drop the fixed planar embedding (Section [6]). 
We first consider biconnected graphs in Section 6.1 and compute cost functions on the interval [0, 3], 
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which can be shown to be convex on that interval, bottom up in the SPQR-tree. In Section [6.2| we 
extend this result to connected graphs using the BC-tree (see Section 2.2 for a definition). 



2 Preliminaries 

In this section we introduce some notations and preliminaries. 

2.1 FlexDraw 

The original FlexDraw problem asks for a given 4-planar graph G = (V, E) with a function 
flex: E — > No U {oo} assigning a flexibility to every edge whether an orthogonal drawing of 
G exists such that every edge e S E has at most flex(e) bends. Such a drawing is called a 
valid drawing of the FlexDraw instance. The problem OptimalFlexDraw is the optimization 
problem corresponding to the decision problem FlexDraw and is defined as follows. Let G = 
(V,E) be a 4-planar graph together with a cost function cost e : No — > M U {cxd} associated with 
every edge e E E having the interpretation that p bends on the edge e cause cost e (p) cost. Then the 
cost of an orthogonal drawing of G is the total cost summing over all edges. A drawing is optimal if 
it has the minimum cost among all orthogonal drawings of G. The task of the optimization problem 
OptimalFlexDraw is to find an optimal drawing of G. 

Since OptimalFlexDraw contains the ATP-hard problem FlexDraw, it is MV-hard itself. 
However, FlexDraw is efficiently solvable for instances with positive flexibility, that is instances in 
which the flexibility of every edge is at least 1. To obtain a similar result for OptimalFlexDraw 
we have to restrict the possible cost functions slightly. 

For a cost function cost e (-) we define the difference function Acost e (-) to be Acost e (p) = 
cost e (p + 1) — cost e (p). A cost function is monotone if its difference function is greater or equal 
to 0. We say that the base cost of the edge e with monotone cost function is b e = cost e (0). The 
flexibility of an edge e with monotone cost function is defined to be the largest possible number 
of bends p for which cost e (p) = b e . As before, we say that an instance G of OptimalFlexDraw 
has positive flexibility if all cost functions are monotone and the flexibility of every edge is positive. 
Unfortunately, we have to restrict the cost functions further to be able to solve OptimalFlexDraw 
efficiently. The cost function cost e (-) is convex, if its difference function is monotone. We call an 
instance of OptimalFlexDraw convex, if every edge has positive flexibility and each cost function 
is convex. Note that this includes that the cost functions are monotone. We provide an efficient 
algorithm solving OptimalFlexDraw for convex instances. 

2.2 Connectivity, BC-Tree and SPQR-Tree 

A graph is connected if there exists a path between any pair of vertices. A separating k-set is a set 
of k vertices whose removal disconnects the graph. Separating 1-sets and 2-sets are cutvertices and 
separation pairs, respectively. A connected graph is biconnected if it does not have a cut vertex 
and triconnected if it does not have a separation pair. The maximal biconnected components of a 
graph are called blocks. The cut components with respect to a separation k-set S are the maximal 
subgraphs that are not disconnected by removing S. 

The block- cutvertex tree (BC-tree) B of a connected graph is a tree whose nodes are the blocks 
and cutvertices of the graph, called B-nodes and C-nodes, respectively. In the BC-tree a block B 
and a cutvertex v are joined by an edge if v belongs to B. If an embedding is chosen for each block, 
these embeddings can be combined to an embedding of the whole graph if and only if B can be 
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Figure 2: The unrooted SPQR-tree of a biconnected planar graph. The nodes fi±, ^ and /X5 are 
P-nodes, \ii is an R-node and \i± is an S-node. The Q-nodes are not shown explicitely. 

rooted at a B-node such that the parent of every other block B in B, which is a cutvertex, lies on 
the outer face of B. 

We use the SPQR-tree introduced by Di Battista and Tamassia [3 [B] to represent all planar 
embeddings of a biconnected planar graph G. The SPQR-tree T of G is a decomposition of G into 
its triconnected components along its split pairs where a split pair is either a separation pair or an 
edge. We first define the SPQR-tree to be unrooted, representing embeddings on the sphere, that 
is planar embeddings without a designated outer face. Let {s,t} be a split pair and let H\ and H2 
be two subgraphs of G such that H\ U H2 = G and H% n H2 = {s, t}. Consider the tree containing 
the two nodes /zi and ^2 associated with the graphs Hi + {s, t} and H2 + {s, t}, respectively. These 
graphs are called skeletons of the nodes fii, denoted by skel(/ij) and the special edge {s,t} is said 
to be a virtual edge. The two nodes /Ui and \i2 are connected by an edge, or more precisely, the 
occurrence of the virtual edges {s,t} in both skeletons are linked by this edge. Now a combinatorial 
embedding of G uniquely induces a combinatorial embedding of skel(//i) and skel(/X2). Furthermore, 
arbitrary and independently chosen embeddings for the two skeletons determine an embedding of 
G, thus the resulting tree can be used to represent all embeddings of G by the combination of all 
embeddings of two smaller planar graphs. This replacement can of course be applied iteratively 
to the skeletons yielding a tree with more nodes but smaller skeletons associated with the nodes. 
Applying this kind of decomposition in a systematic way yields the SPQR-tree as introduced by 
Di Battista and Tamassia [3 |5]. The SPQR-tree T of a biconnected planar graph G contains four 
types of nodes. First, the P-nodes having a bundle of at least three parallel edges as skeleton 
and a combinatorial embedding is given by any ordering of these edges. Second, the skeleton of 
an R-node is triconnected, thus having exactly two embeddings [18] . and third, S- nodes have a 
simple cycle as skeleton without any choice for the embedding. Finally, every edge in a skeleton 
representing only a single edge in the original graph G is formally also considered to be a virtual 
edge linked to a Q-node in T representing this single edge. Note that all leaves of the SPQR-tree T 
are Q-nodes. Besides from being a nice way to represent all embeddings of a biconnected planar 
graph, the SPQR-tree has only size linear in G and Gutwenger and Mutzel [TU] showed how to 
compute it in linear time. Figure [2] shows a biconnected planar graph together with its SPQR-tree. 

Often the SPQR-tree T of a biconnected planar graph G is assumed to be rooted in a Q-node 
representing all planar embeddings with the corresponding edge on the outer face. In contrast to 
previous results, we assume the SPQR-tree T to be rooted in some node r, which may be a Q-node 
or an inner node. In the following we describe the interpretation of the SPQR-tree with root r. 
Every node fi, apart form r itself, has a unique parent and thus its skeleton skel(/i) contains a 
virtual edge corresponding to this parent. We refer to this virtual edge as the parent edge. A 
planar embedding £ of G is represented by T with root r if the embedding induced on the skeleton 
skel(/i) of every node /i/r has the parent edge on the outer face. The embedding of skel(r) is not 
restricted, thus the choice of the outer face makes a difference for the root. 
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For every node \i in the SPQR-tree T apart from the root r we define the pertinent graph of /j,, 
denoted by pert(^), as follows. The pertinent graph of a Q-node is the edge associated to it. The 
pertinent graph of an inner node [i is recursively defined to be the graph obtained by replacing 
all virtual edges apart from the parent edge by the pertinent graphs of the corresponding children 
in T. The expansion graph of a virtual edge e in skel(/u) is the pertinent graph of fjf where fjf is 
the child of \i corresponding to the virtual edge e with respect to the root /j,. 

2.3 Orthogonal Representation 

Two orthogonal drawings of a 4-planar graph G are equivalent, if they have the same topology, that 
is the same planar embedding, and the same shape in the sense that the sequence of right and left 
turns is the same in both drawings when traversing the faces of G. To make this precis, we define 
orthogonal representations, originally introduced by Tamassia [15], as equivalence classes of this 
equivalence relation between orthogonal drawings. To ease the notation we first only consider the 
biconnected case. 

Let r be an orthogonal drawing of a biconnected 4-planar graph G. In the planar embedding £ 
induced by Y every edge e is incident to two different faces, let / be one of them. When traversing 
/ in clockwise order (counter-clockwise if / is the outer face) e may have some bends to the right 
and some bends to the left. We define the rotation of e in the face / to be the number of bends to 
the right minus the number of bends to the left and denote the resulting value by rot(e/). Similarly, 
every vertex v is incident to several faces, let / be one of them. Then we define the rotation of v 
in /, denoted by rot(uj), to be 1, — 1 and if there is a turn to the right, a turn to the left and 
no turn, respectively, when traversing / in clockwise direction (counter-clockwise if / is the outer 
face). The orthogonal representation 1Z belonging to T consists of the planar embedding £ of G 
and all rotation values of edges and vertices, respectively. It is easy to see that every orthogonal 
representation has the following properties. 

(I) For every edge e incident to the faces f± and /2 the equation rot(e/ 1 ) = — rot(e/ 2 ) holds. 
(II) The sum over all rotations in a face is 4 for inner faces and —4 for the outer face. 
(Ill) The sum of rotations around a vertex v is 2 • (deg(i>) — 2). 

Tamassia showed that the converse is also true [T5], that is 7Z is an orthogonal representation 
representing a class of orthogonal drawings if the rotation values satisfy the above properties. He 
moreover describes a flow network such that every flow in the flow network corresponds to an 
orthogonal representation. A modification of this flow network can also be used to solve Op- 
timalFlexDraw but only for the case that the planar embedding is fixed. In some cases we 
also write Tot-ji(-) instead of rot(-) to make clear to which orthogonal representation we refer to. 
Moreover, the face in the index is sometimes omitted if it is clear which face is meant. 

When extending the term orthogonal representation to not necessarily biconnected graphs there 
are two differences. First, a vertex v with deg(t;) = 1 may exist. Then v is incident to a single face 
/ and we define the rotation rot(uj) to be —2. Note that the rotations around every vertex v still 
sum up to 2 • (deg(w) — 2) . The second difference is that the notation introduced above is ambiguous 
since edges and vertices may occur several times in the boundary of the same face. For example 
a bridge e is incident to the face / twice, thus it is not clear which rotation is meant by rot(e/). 
However, it will always be clear from the context, which incidence to the face / is meant by the 
index /. Thus, we use for connected graphs the same notation as for biconnected graphs. 

Let G be a 4-planar graph with orthogonal representation 1Z and two vertices s and t incident 
to a common face /. We define vrj(s, t) to be the unique shortest path from s to t on the boundary 
of /, when traversing / in clockwise direction (counter-clockwise if / is the outer face). Let 



6 



-o — o 



-o — o 



-o — o 



o — o 6 — 6 — o 6 — 6 



O — 6 — 6 — Q — 6 



o — o 



-6 — o 



? 



Figure 3: On the left three tight orthogonal drawings are stacked together. This is not possible on 
the right side, since the black vertices have angles larger than 90° in internal faces. 



v±, . . . , Vk = t be the vertices on the path nf(s, t). The rotation of n(s, t) is defined as 

k-1 k-1 

rot(7r(s, t)) = rot ({^> v i+i}) + rot ( 7 



[Vi 



where all rotations are with respect to the face /. 

Note that it does not depend on the particular drawing of a graph G how many bends each edge 
has but only on the orthogonal representation. Thus we can continue searching for valid and optimal 
orthogonal representations instead of drawings to solve FlexDraw and OptimalFlexDraw, 
respectively. 

Let G be a 4-planar graph with positive flexibility and valid orthogonal representation 1Z and 
let {s,t} be a split pair. Let further if be a split component with respect to {s,t} such that the 
orthogonal representation S of H induced by TZ has {s,t} on the outer face /. The orthogonal 
representation S of H is called tight with respect to the vertices s and t if the rotations of s and 
t in internal faces are 1, that is s and t form 90°-angles in internal faces of H. Blasius et al. [2J 
Lemma 2] show that S can be made tight with respect to s and t, that is there exists a valid tight 
orthogonal representation of H that is tight. Moreover, this tight orthogonal representation can 
be plugged back into the orthogonal representation of the whole graph G. We call an orthogonal 
representation 1Z of the whole graph G tight, if every split component having the corresponding 
split pair on its outer face is tight with respect to its split pair. It follows that we can assume 
without loss of generality that every valid orthogonal representation is tight. This has two major 
advantages. First, if we have for example a chain of graphs and orthogonal representations of 
each graph in the chain, we can combine these orthogonal representations by simply stacking them 
together; see Figure [3j Note that this may not be possible if the orthogonal representations are not 
tight. Second, the shape of the outer face / of a split component with split pair {s, t} is completely 
determined by the rotation of 717 (s, t) and the degrees of s and t, since the rotation at the vertices s 
and t in the outer face only depends on their degrees. In the following we assume every orthogonal 
representation to be tight. 



2.4 Flow Network 

A cost flow network (or flow network for short) is a tuple N = (V, A, COST, dem) where (V, A) is 
a directed (multi-)graph, COST is a set containing a cost function cost a : No — > KU {00} for each 
arc a £ A and dem : V — > Z is the demand of the vertices. A flow in N is a function <fi : A — > No 
assigning a certain amount of flow to each arc. A flow eft is feasible, if the difference of incoming 
and outgoing flow at each vertex equals its demand, that is 

dem('w) = 4>{u, v) — <p(v, u) for all v E V. 

(u,v)EA (v,u)eA 

The cost of a given flow (j) is the total cost of the arcs caused by the flow 4>, that is 

cost(0) = cost a (0(a)). 

aeA 
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A feasible flow <p in N is called optimal if cost(^) < cost(^>') holds for every feasible flow (j)'. 

If the cost function of an arc a is on an interval [0, c] and oo on (c, oo), we say that a has 
capacity c. 

A flow network TV" is called convex if the cost functions on its arcs are convex. In the flow 
networks we consider, every arc a E A has a corresponding arc a' £ A between the same vertices 
pointing in the opposite direction. A flow (p is normalized if (j)(a) = or (ft(a') = for each of these 
pairs. Since we only consider convex flow networks a normalized optimal flow does always exist. 
Thus we assume without loss of generality that all flows are normalized. We simplify the notation 
as follows. If we talk about an amount of flow on the arc a that is negative, we instead mean the 
same positive amount of flow on the opposite arc a' . In many cases minimum-cost flow networks are 
only considered for linear cost functions, that is each unit of flow on an arc causes a constant cost 
defined for that arc. Note that the cost functions in a convex flow network iV are piecewise linear 
and convex according to our definition. Thus, it can be easily formulated as a flow network with 
linear costs by splitting every arc into multiple arcs, each having linear costs. It is well known that 
flow networks of this kind can be solved in polynomial time. The best known running time depends 
on additional properties that N may satisfy. We use an algorithm computing a minimum-cost flow 



in the network N as black box and denote the necessary running time by Ifl ow (|iV|). In Section 6.3 
we have a closer look on which algorithm to use. 

Let u, v S V be two nodes of the convex flow network N with demands dem(«) and dem(u). 
The parameterized flow network with respect to the nodes u and v is defined the same as N but 
with a parameterized demand of dem(u) — p for u and dem(w) + p for v where p is a parameter. The 
cost function costAr(p) of the parameterized flow network N is defined to be cost(^) of an optimal 
flow (ft in N with respect to the parameterized demands determined by p. Note that increasing 
p by 1 can be seen as pushing one unit of flow from u to v. We define the optimal parameter po 
to be the parameter for which the cost function is minimal among all possible parameters. The 
correctness of the minimum weight path augmentation method to compute flows with minimum 
costs implies the following theorem [7|. 

Theorem 1. The cost function of a parameterized flow network is convex on the interval [po,oo], 
where po is the optimal parameter. 

Proof. Let N = (V, A, COST, dem) be a parameterized flow network and let (fio be a minimum-cost 
flow in N with respect to the optimal parameter p$. To simplify notation, we assume po = 0. The 
residual network Rq with respect to 4>o is the graph (V, A) with a constant cost costo(a) assigned to 
every arc a such that costo(a) is the amount of cost in N that has to be payed to push an additional 
unit of flow along a, with respect to the given flow <f)Q. Note that this cost may be negative. It is 
well known that an optimal flow <j>i with respect to the parameter 1 can be computed by pushing 
one unit of flow along a path from u to v with minimum weight in Rq [7i. Moreover, we can 
continue and compute an optimal flow (f>k+i by augmenting <pk along a minimum weight path in 
the residual network with respect to the flow 0^. Assume we augment 0& along the path ir^ 
causing cost costfc(7Tfc) to obtain an optimal flow 4>k+i with respect to the parameter k + 1 and 
then we augment along a path Tik+i in Rk+i with cost costfc +1 (7Tfc +1 ) to obtain an optimal flow 
4>k+2 with respect to the parameter k + 2. To obtain the claimed convexity we have to show that 
cost fc (7r fc ) < cost /c+1 (7r A . +1 ) holds. 

If 7Tfc and Tik+i contain an arc a in the same direction, then costfc(a) < costfc+i(a) holds by the 
convexity of the cost function of a. If irk contains the arc a and n^+i contains the arc a' in the 
opposite direction then costfc(a) = — costfc+i(a') holds. Assume irk and nt+i share such an arc in 
the opposite direction. Then we remove this arc in both directions, splitting each of the paths 7r& 
and 7Tfc +1 into two subpaths. We define two new paths tt and n' by concatenating the first part of ir^ 
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Figure 4: Since a strictly directed path from t to s has a lower bound for its rotation this yields 
upper bounds for paths from s to t (Lemma [TJ) . 

with the second part of 7Tk+i and vice versa, respectively. This can be done iteratively, thus we can 
assume that tt and tt' do not share arcs in the opposite direction. We consider the cost of tt and tt' in 
the residual network Rk ■ Obviously, for an arc a that is exclusively contained either in tt or in tt' we 
have costfc(a) = costfc +1 (a). For an arc that is contained in tt and tt' we have costfc(a) < costfc+i(a). 
Moreover, for every pair of arcs a and a' that was removed we have costfc(a) = — costfc+i(a'). This 
yields the inequality costfc(-7Tfc) + cost/ c+ i(7Tj! c+1 ) > costfc(7r) + costfc(7r'). Since tt^ was a path with 
smallest possible weight in Rk we have costfc(vrfc) < costfe(7r) and cost&(7Tfc) < costfc(vr'). With the 
above inequality this yields costfc + i(7Tfc + i) > costfc(7Tfc). □ 



3 Valid Drawings with Fixed Planar Embedding 

In this section we consider the problem FlexDraw for the case that the planar embedding is fixed. 
We show that the existence of a valid orthogonal representation implies the existence of a valid 
orthogonal representation with special properties. We first show the following. Given a biconnected 
4-planar graph with positive flexibility and an orthogonal representation 1Z such that two vertices s 
and t lie on the outer face /, then the rotation along irt (s, t) can be reduced by 1 if it is at least 0. 
This result is a key observation for the algorithm solving the decision problem FlexDraw [2|. It 
in a sense shows that "rigid" graphs that have to bent strongly do not exists. This kind of graphs 
play an important role in the AA-P-hardness proof of O-embeddability by Garg and Tamassia [3]. 
Moreover, we show the existence of a valid orthogonal representation 1Z' inducing the same planar 
embedding and having the same angles around vertices as 1Z such that every edge has at most three 
bends in 1Z' , except for a single edge on the outer face with up to five bends. If we allow to change 
the embedding slightly, this special edge has only up to four bends. 

Let G be a 4-planar graph with positive flexibility and valid orthogonal representation 1Z, and 
let e be an edge. If the number of bends of e equals its flexibility, we orient e such that its bends 
are right bends. Otherwise, e remains undirected. We define a path tt = [v\, . . . ,Vk) in G to be a 
directed path, if the edge {vi, v; L+ i} (for i G {1, . . . , k — 1}) is either undirected or directed from Vi 
to Uj-fi. A path containing only undirected edges can be seen as directed path for both possible 
directions. The path tt is strictly directed, if it is directed and does not contain undirected edges. 
These terms directly extend to (strictly) directed cycles. Given a (strictly) directed cycle C the 
terms left(C) and right(C) denote the set of edges and vertices of G lying to the left and right of C, 
respectively, with respect to the orientation of C. A cut (U, V\U) is said to be directed from U to 
V \ U, if every edge {u, v} with u G U and v £ V \ U is either directed from u to v or undirected. 
According to the above definitions a cut is strictly directed from U to V \ U if it is directed and 
contains no undirected edges. Before we show how to unwind an orthogonal representation that is 
bent strongly we need the following technical lemma. 

Lemma 1. Let G be a graph with positive flexibility and vertices s and t such that G + st is 
biconnected and J^-planar. Let further 1Z be a valid orthogonal representation with s and t incident 
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to the common face f such that TTf(t, s) is strictly directed from t to s. Then the following holds. 

(1) iotn(iTf(s, t)) < —3 if f is the outer face and G does not consist of a single path 

(2) votTi(wf(s,t)) < — 1 if f is the outer face 

(3) vot n (TTf(s,t)) <5 

Proof We first consider the case where / is the outer face (Figure Qa)), that is cases (1) and (2). 
Due to the fact that TTf(t, s) is strictly directed from t to s and the flexibility of every edge is positive, 
each edge on TTf(t,s) has rotation at least 1. Moreover, the rotations at vertices along the path 
TTf(t, s) are at least —1 since 7r/(t, s) is simple as G + st is biconnected. Since the number of internal 
vertices on a path is one less than the number of edges this yields Tot(iVf(t, s)) > 1; see Figure |4^b) . 
If G consists of a single path this directly yields rot(7i/(,s, t)) < —1 and thus concludes case (2). For 
case (1) first assume that the degrees of s and t are not 1 (Figure |4|b) ) , that is rot(sj), rot(i/) € 
{— 1, 0, 1} holds. Since / is the outer face the equation rot(7T/(s, i))+rot(iy)+rot(7ry(i, s))+rot(s/) = 
—4 holds and directly implies the desired inequality rot(-7r/(s, t)) < —3. In the case that for 
example t has degree 1 (and deg(s) > 0), we have rot(tj) = —2 and rot(sj) 6 { — 1,0, 1}, thus the 
considerations above only yield rot(7rj(s, t)) < —2. However, in this case there necessarily exists a 
vertex t' where the paths 7Tf(s,t) and 7Tf(t,s) split, as illustrated in Figure |4^c). More precisely, 
let t' be the first vertex on nf(s,t) that also belongs to 7T/(i,s). Obviously, the degree of t' is 
at least 3 and thus rot(ij) (with respect to the path irf(t,s)) is at least 0. Hence we obtain the 
stronger inequality vot(n f(t, s)) > 2 yielding the desired inequality rot(7ry(s, t)) < —3. If s and t 
both have degree 1 we cannot only find the vertex t! but also the vertex s' where the paths tt/(s, t) 
and 7Tf(t,s) split. Since G + st is biconnected these two vertices are distinct and the estimation 
above works, finally yielding rot(7T/(s, t)) < —3. 

If / is an internal face (Figure |4^d)), that is case (3) applies, we start with the equation 
rot(7T/-(s, t)) + rot(t/) + rot(7T/(t, s)) + rot(s/) = 4. First we consider the case that neither t nor s 
have degree 1. Thus, rot(fy), rot(sy) E {—1,0,1}. With the same argument as above we obtain 
rot(7T/(i, s)) > 1 and hence rot(-7r/(s, t)) < 5; see Figure |4](e) . Now assume that t has degree 1 and s 
has larger degree. Then rot(tj) = —2 holds and the above estimation does not work anymore. 
Again, at some vertex t' the paths ?r/(i, s) and vrj(s, t) split as illustrated in Figure |4|T) . Obviously, 
the degree of t' needs to be greater than 2 and thus rot(tj) is at least 0. This yields iot(irf(t, s)) > 2 
in the case that deg(t) = 1, compensating rot(i/) = —2 (instead of rot(ij) > — 1 in the other case). 
To sum up, we obtain the desired inequality rot(7rj (s, t)) < 5. The case deg(s) = deg(t) = 1 works 
analogously. □ 

The flex graph G^ of G with respect to a valid orthogonal representation 7Z is defined to be the 
dual graph of G such that the dual edge e* is undirected if e is undirected, otherwise it is directed 
from the face right of e to the face left of e. Figure [5^a) shows an example graph with an orthogonal 
drawing together with the corresponding flex graph. Assume we have a simple directed cycle C 
in the flex graph. Then bending along this cycle yields a new valid orthogonal representation TV 
which is defined as follows. Let e* = (/i, $2) be an edge contained in C dual to e. Then we decrease 
rot(ej 1 ) and increase rot(ej 2 ) by 1. It can be easily seen that the necessary properties for 1Z' to be 
an orthogonal representation are satisfied. Obviously, rot7^./(e/ 1 ) = — rot-R/(e/ 2 ) holds and rotations 
at vertices did not change. Moreover, the rotation around a face / does not change since / is either 
not contained in C or it is contained in C, but then it has exactly one incoming and exactly one 
outgoing edge. Note that bending along a cycle in the flex graph preserves the planar embedding 
of G and for every vertex the rotations in all incident faces. The following lemma shows that a 
high rotation along a path irf(s,t) for two vertices s and t sharing the face / can be reduced by 1 
using a directed cycle in the flex graph. 
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Figure 5: (a) An orthogonal representation and the corresponding flex graph where every edge has 
flexibility 1. (b, c, d) Illustration of Lemma [2] 

Lemma 2. Let G be a biconnected J^-planar graph with positive flexibility, a valid orthogonal rep- 
resentation 1Z and s and t on a common face f. The flex graph G^ contains a directed cycle C 
such that f £ C, s £ left(C) and t £ right(C), if one of the following conditions holds. 

(1) votTi(7r f(s,t)) > —2, / is the outer face and TTf(s,t) is not strictly directed from t to s 

(2) rot-ji(TTf(s,t)) > and f is the outer face 

(3) rot 7 j( 7 r / ( S ,t))>6 

Proof. Figure |5|b) shows the path 7T/(s,i) together with the desired cycle C. Due to the duality 
of a cycle in the dual and a cut in the primal graph a directed cycle C in G^ having s and t to 
the left and to the right of C, respectively, induces a directed cut in G that is directed from s to 
t and vice versa. Recall that directed cycles and cuts may also contain undirected edges. Assume 
for contradiction that such a cycle C does not exist. 

Claim 1. The graph G contains a strictly directed path tt from t to s. 

Every cut (S, T) with T = V \ S, s £ S and t £ T separating s from t must contain an edge that is 
directed from T to S, otherwise this cut would correspond to a cycle C in the flex graph that does 
not exist by assumption. Let T be the set of vertices in G that can be reached by strictly directed 
paths from t. If T contains s we found the path tt strictly directed from t to s. Otherwise, (S, T) 
with S = V \ T is a cut separating S from T and there cannot be an edge that is directed from a 
vertex in T to a vertex in S which is a contradiction, and thus the path tt strictly directed from t 
to s exists, which concludes the proof of the claim. 

Let G' be the subgraph of G induced by the paths tt and TTf(s,t) together with the orthogonal 
representation 1Z' induced by 1Z. 

We first consider case (1). Let /' be the outer face of the orthogonal representation 7Z'. Obvi- 
ously, TTf/(s,t) = TTf(s,t) and tt = TTf(t,s) holds, see Figure |5|c). Moreover, the graph G' + st is 
biconnected and G' does not consist of a single path since TTf(s, t) and TTf(t, s) are different due to 
the assumption that Wf(s,t) is not strictly directed from t to s. Since TTf(t,s) is strictly directed 
from t to s we can use Lemma [T|l) yielding rot-R/(7i//(s, t)) < —3 and thus rot-ft(7T/(s, t)) < —3, 
which is a contradiction. 

For case (2) exactly the same argument holds except for the case where the strictly directed 
path tt is the path 7Tf(s,t) strictly directed from t to s. In this case we have to use Lemma [T|2) 
instead of Lemma [T]^l) yielding Tot-ji(TTf(s,t)) < —1, which is again a contradiction. 

In case (3) the subgraph G' of G induced by the two paths tt and 7Ty(s, t) again contains s and t 
on a common face /', which may be the outer or an inner face, see Figure [5jc) and Figure [5jd), 
respectively. In both cases we obtain rotTi(TTf(s, t)) < 5 due to Lemma[T]^3), which is a contradiction. 

□ 

Lemma [2] directly yields the following corollary, showing that graphs with positive flexibility 
behave very similar to single edges with positive flexibility. 
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Corollary 1. Let G be a graph with positive flexibility and vertices s and t such that G + st is 
biconnected and ^-planar. Let further 1Z be a valid orthogonal representation with s and t on the 
outer face f such that p = rotn(ir f(s,t)) > 0. For every rotation p' £ [—l,p] there exists a valid 
orthogonal representation 1Z' with Tot-ji'(7rf(s,t)) = p' . 

Proof. For the case that G itself is biconnected, the claim follows directly from Lemma|2]^2) , since we 
can reduce the rotation along 7i/(s, t) stepwise by 1, starting with the orthogonal representation 1Z, 
until we reach a rotation of —1. For the case that G itself is not biconnected we add the edge {s,t} 
to the orthogonal representation 1Z such that the path 7r/(s,i) does not change, that is 7Tf(t,s) 
consists of the new edge Again Lemma [2^2) can be used to reduce the rotation stepwise 

down to —1. □ 

As edges with many bends imply the existence of paths with high rotation, we can use Lemma[2] 
to successively reduce the number of bends of every edge down to three, except for a single edge 
on the outer face. Since we only bend along cycles in the flex graph, neither the embedding nor 
the angles around vertices are changed. 

Theorem 2. Let G be a biconnected J^-planar graph with positive flexibility, having a valid orthogo- 
nal representation. Then G has a valid orthogonal representation with the same planar embedding, 
the same angles around vertices and at most three bends per edge, except for at most one edge on 
the outer face with up to five bends. 

Proof. In the following we essentially pick an edge with more than three bends, reduce the number 
of bends by one and continue with the next edge. After each of these reduction steps we set the 
flexibility of every edge down to max{/>, 1}, where p is the number of bends it currently has. This 
ensures that in the next step the number of bends of each edge either is decreased, remains as it is 
or is increased from zero to one. 

We start with an edge e = {s,t} that is incident to two faces f\ and f% and has more than 
three bends. Due to the fact that we traverse inner faces in clockwise and the outer face in counter- 
clockwise direction, the edge e forms in one of the two faces the path from s to t and in the other face 
the path from t to s. Assume without loss of generality that 71"^ (t, s) and vrj 2 (s, t) are the paths on 
the boundary of fi and /2, respectively, that consist of e. Note that rot(-7rj 1 (t, s)) = — rot(7r/ 2 (s, i)) 
holds and we assume that rot(7Tj 1 (t, s)) is not positive. As e was assumed to have more than three 
bends, the inequality rot(7r/ 1 (t, s)) < —4 holds. We distinguish between the two cases that f\ is an 
inner or the outer face. We first consider the case that f\ is an inner face; Figure [6ja) illustrates 
this situation for the case where e has four bends. Then the rotations around the face f\ sum up 
to 4. As the rotations at the vertices s and t can be at most 1, we obtain rot(7Tj 1 (s, t)) > 6. Thus 
we can apply Lemma [2^3) to reduce the rotation of irf 1 (s,t) by bending along a cycle in the flex 
graph that contains f\ and separates s from t. Obviously, this increases the rotation along TTf^t, s) 
by 1 and thus reduces the number of bends of e by 1. 

For the case that f\ is the outer face we first ignore the case where e has four or five bends and 
show how to reduce the number of bends to five; Figure [6^b) shows the case where e has six bends. 
Thus the inequality rot(-7r/ 1 (i, s)) < —6 holds. As the rotations around the outer face f\ sum up 
to —4 and the rotations at the vertices s and t are at most 1, the rotation along ^(s,^ must be 
at least 0. Thus we can apply Lemma [2^2) to reduce the rotation of ^^(sjt) by 1, increasing the 
rotation along 71^ (i, s), and thus reducing the number of bends of e by one. 

Finally, we obtain an orthogonal representation having at most three bends per edge except for 
some edges on the outer face with four or five bends having their negative rotation in the outer face. 
If there is only one of these edges left we are done. Otherwise let e = {s, t} be one of the edges with 
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Figure 6: Reducing the number of bends on edges (Theorem 2 



rot(-7Tj(t, s)) £ {—5, —4}, where / is the outer face. Then the inequality rot(7r/(s, t)) > —2 holds by 
the same argument as before and we can apply Lemma [2^1) to reduce the rotation, if we can ensure 
that 7Tf(s, t) is not strictly directed from t to s. To show that, we make use of the fact that tt/(s, t) 
contains an edge e' = {u, v} with at least four bends due to the assumption that e was not the 
only edge with more than three bends. Assume without loss of generality that u occurs before v 
on TTf(s,t), thus iTf(s,t) splits into the three parts iTf(s,u), irf(u,v) and iTf(v,t). Recall that 
rot(7Tj(s, t)) > —2 holds and thus rot(7Tj(s, u)) +rot(u) +rot(7T/(u, v)) +rot(t> ) + Tot(iTf(v, t)) > —2. 
As the rotation at the vertices u and v is at most 1 and the rotation of iTf{u,v) at most —4 it 
follows that rot(7rj(s, u)) + rot(7rj(u, t)) > 0. Figure [6]^c) illustrates the situation for the case 
where e and e' have four bends and rot(7Tj(s, u)) = rot(-K j (v , t)) = 0. Note that at least one of the 
two paths is not degenerate in the sense that s ^ u or v ^ t, otherwise the total rotation around the 
outer face would be at most —6, which is a contradiction. Assume without loss of generality that 
rot(-7rj(s, u)) > 0. It follows that nj(s, u) cannot be strictly directed from u to s and since 7i/(s, u) 
is a subpath of irf(s,t) the path 7T/(s,t) cannot be strictly directed from t to s. This finally shows 
that we can use part (1) of Lemma [2] implying that we can find a valid orthogonal representation 
such that at most a single edge with four or five bends remains, whereas all other edges have at 
most three bends. □ 

If we allow the embedding to be changed slightly, we obtain an even stronger result. Assume 
the edge e lying on the outer face has more than three bends. If e has five bends, we can reroute 
it in the opposite direction around the rest of the graph, that is we can choose the internal face 
incident to e to be the new outer face. In the resulting drawing e has obviously only three bends. 
Thus the following result directly follows from Theorem [2} 

Corollary 2. Let G be a biconnected 4-planar graph with positive flexibility having a valid orthogo- 
nal representation. Then G has a valid orthogonal representation with at most three bends per edge 
except for possibly a single edge on the outer face with four bends. 

Note that Corollary [2] is restricted to biconnected graphs. For general graphs it implies that 
each block contains at most a single edge with up to four bends. Figure [7] illustrates an instance 
of FlexDraw with linearly many blocks and linearly many edges that are required to have four 
bends, showing that Corollary [2] is tight. 

Theorem [2] implies that it is sufficient to consider the flexibility of every edge to be at most 5, 
or in terms of costs we want to optimize, it is sufficient to store the cost function of an edge only 
in the interval [0, 5]. However, there are two reasons why we need a stronger result. First, we want 
to compute cost functions of split components and thus we have to limit the number of "bends" 
they can have (see the next section for a precise definition of bends for split components). Second, 
as mentioned in the introduction (see Figure [TJ the cost function of a split component may already 
be non-convex on the interval [0,5]. Fortunately, the second reason is not really a problem since 
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Figure 7: An instance of FlexDraw requireing linearly many edges to have four bends. Flexibilites 
are 1 except for the thick edges with flexibility 4. 

there may be at most a single edge with up to five bends, all remaining edges have at most three 
bends and thus we only need to consider their cost functions on the interval [0, 3] . 

In the following section we focus on dealing with the first problem and strengthen the results so 
far presented by extending the limitation on the number of bends to split components. Note that 
a split pair inside an inner face of G with a split component H having a rotation less than —3 on 
its outer face implies a rotation of at least 6 in some inner face of G. Thus, we can again apply 
Lemma (2^3) to reduce the rotation showing that split components and single edges can be handled 
similarly. However, by reducing the rotation for one split component, we cannot avoid that the 
rotation of some other split component is increased. For single edges we did that by reducing the 
flexibility to the current number of bends. In the following section we extend this technique by 
defining a flexibility not only for edges but also for split components. We essentially show that all 
results we presented so far still apply, if we allow this kind of extended flexibilities. 

4 Flexibility of Split Components and Nice Drawings 

Let G be a biconnected 4-planar graph with SPQR-tree T and let T be rooted at some node r. 
Recall that we do not require r to be a Q-node. Let \x be a node of T that is not the root r. Then fj, 
has a unique parent and skel(//) contains a unique virtual edge e = {s,t} that is associated with 
this parent. We call the split-pair {s,t} a principal split pair and the pertinent graph pert(/i) with 
respect to the chosen root a principal split component. The vertices s and t are the poles of this 
split component. Note that a single edge is also a principal split component except for the case 
that its Q-node is chosen to be the root. A planar embedding of G is represented by T with the 
root r if the embedding of each skeleton has the edge associated with the parent on the outer face. 

Let 1Z be a valid orthogonal representation of G such that the planar embedding of 1Z is 
represented by T rooted at r. Consider a principal split component H with respect to the split 
pair {s, t} and let S be the orthogonal representation of H induced by 1Z. Note that the poles 
s and t are on the outer face / of S. We define max{| iots(irf(s, t))\, \ icots(7ff(t, s))\} to be the 
number of bends of the split component H. Note that this is a straightforward extension of the 
term bends as it is used for edges. With this terminology we can assign a flexibility flex(iif) to 
a principal split component H and we define the orthogonal representation 1Z of G to be valid if 
and only if H has at most iiex(H) bends. We say that the graph G has positive flexibility if the 
flexibility of every principal split component is at least 1 , which is straightforward extension of the 
original notion. 

We define a valid orthogonal representation of G to be nice if it is tight and if there is a root r 
of the SPQR-tree such that every principal split component has at most three bends and the edge 
corresponding to r in the case that r is a Q-node has at most five bends. The main result of this 
section will be the following theorem, which directly extends Theorem [2j 

Theorem 3. Every biconnected J^-planar graph with positive flexibility having a valid orthogonal 
representation has an orthogonal representation with the same planar embedding and the same 
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Figure 8: Augmentation of G by the safety edges en(s,t) and eu(t,s). 



angles around vertices that is nice with respect to at least one node chosen as root of its SPQR-tree. 

Before we prove Theorem [3] we need to make some additional considerations. In particular we 
need to extend the flex-graph such that it takes the flexibilities of principal split components into 
account. The extended version of the flex graph can then be used to obtain a result similar to 
Lemma [2] which was the main tool to proof Theorem [2j Another difficulty is that it depends on 
the chosen root which split components are principal split components. For the moment we avoid 
this problem by choosing an arbitrary Q-node to be the root of the SPQR-tree T ■ Thus we only 
have to care about the flexibilities of the principal split components with respect to the chosen 
root. One might hope that the considerations we make for the flex- graph in the case of a fixed root 
still work, if we consider the principal split components with respect to all possible roots at the 
same time. However, this fails as we will see later, making it necessary to consider internal vertices 
as the root. 

Assume that the SPQR-tree T of G is rooted at the Q-node corresponding to an arbitrary 
chosen edge. Let H be a principal split component with respect to the chosen root with the poles s 
and t. In the embedding of G the outer face / of H splits into two faces f\ and fa, where the 
path -Kf(s,t) is assumed to lie in f\ and nf(t,s) is assumed to lie in /a, that is ^^(s^t) = nf(s,t) 
and TTf 2 (t, s) = nf(t, s). We augment G by inserting the edge {s, t} twice, embedding one of them 
in fx and the other in fi- We denote the edge {s, t} inserted into the face f\ by en(s,t) and the 
edge inserted into /2 by e#(i, s). Figure [8] illustrates this process and shows how the dual graph 
of G changes. We call the new edges ejj(s, t) and eu (t, s) safety edges and define the extended flex 
graph G x as before, ignoring that some edges have a special meaning. To simplify notation we often 
use the term flex graph, although we refer to the extended flex graph. Note that every cycle in the 
flex graph that separates s from t and thus crosses ir(s, t) and ir(t, s) needs to also cross the safety 
edges e#(s, t) and e#(i, s). Thus we can use the safety edges to ensure that the flex graph respects 
the flexibility of H by orienting them if necessary. More precisely, we orient the safety edge en(s,t) 
from t to s if rot(7r(s, t)) = — flex(ii~) and similarly e#(i, s) from s to t if rot(7r(t, s)) = — flex(iJ). 
This ensures that the rotations along ir(s, t) and ir(t, s) cannot be reduced below — flex(ff) by 
bending along a cycle in the flex graph. Moreover, rot(7r(s,t)) cannot be increased above Hex(H) 
as otherwise rot(7r(i, s)) has to be below — Qex(H) and vice versa. To sum up, we insert the safety 
edges next to the principal split component H and orient them if necessary to ensure that bending 
along a cycle in the flex graph respects not only the flexibilities of single edges but also the flexibility 
of the principal split component H. 

Since adding the safety edges for the graph H is just a technique to respect the flexibility 
of H by bending along a cycle in the flex graph, we do not draw them. Note that the augmented 
graph does not have maximum degree 4 anymore but this is not a problem since we do not draw 
the safety edges. However, we formally assign an orthogonal representation to the safety edges 
by essentially giving them the shape of the paths they "supervise". More precisely, the edges 
en(s,t) and e#(i, s) have the same rotations as the paths 7r(s,t) and ir(t, s) on the outer face of 
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H, respectively. Moreover, the angles at the vertices s and t are also assumed to be the same as 
for these two paths. 

As we do not only want to respect the flexibility of a single split component, we add the safety 
edges for each of the principal split components at the same time. Note that the augmented graph 
remains planar as we only add the safety edges for the principal split components with respect 
to a single root. It follows directly that the considerations above still work, which would fail if 
the augmented graph was non-planar. This is the reason why we cannot consider the principal 
split components with respect to all roots at the same time. The following lemma directly extends 
Lemma [2] to the case where the extended flex graph is considered. 

Lemma 3. Let G be a biconnected J^-planar graph with positive flexibility, a valid orthogonal rep- 
resentation 1Z and s and t on a common face f. The extended flex graph G^ contains a directed 
cycle C such that f £ C , s £ left(C) and t £ right(C), if one of the following conditions holds. 

(1) votTi(wf(s,t)) > —2, f is the outer face and nf(s,t) is not strictly directed from t to s 

(2) votTi(wf(s,t)) > and f is the outer face 

(3) rot^(vr / (s,t)) > 6 

Proof. As in the proof of Lemma [2] we assume for contradiction that the cycle C does not exists, 
yielding a strictly directed path from t to s in G. This directly yields the claim, if we can apply 
Lemma [T] as before. The only difference to the situation before is that the directed path from t 
to s may contain some of the safety edges. However, by definition a safety edge eu(u, v) is directed 
from v to u if and only if Tot(ir(u,v)) = —fiex(H). As flex(iJ) is positive rot(-7r(u, v)) has to be 
negative and thus the rotation along en(u,v) when traversing it from v to u is at least 1. Thus, 
it does not make a difference whether the directed path from t to s consists of normal edges or 
may contain safety edges. Hence, Lemma [T] extends to the augmented graph containing the safety 
edges, which concludes the proof. □ 

Now we are ready to prove Theorem [3j To improve readability we state it again. 

Theorem 3. Every biconnected 4-planar graph with positive flexibility having a valid orthogonal 
representation has an orthogonal representation with the same planar embedding and the same 
angles around vertices that is nice with respect to at least one node chosen as root of its SPQR-tree. 

Proof. Let 1Z be a valid orthogonal representation of G. We assume without loss of generality 
that 1Z is tight. Since the operations we apply to 1Z in the following do not affect the angles 
around vertices, the resulting orthogonal representation is also tight. Thus it remains to enforce 
the more interesting condition for orthogonal representations to be nice, that is reduce the number 
of bends of principal split components down to three. As mentioned before, the SPQR-tree T of 
G is initially rooted at an arbitrary Q-node. Let e re f be the corresponding edge. As in the proof 
of Theorem [2] we start with an arbitrary principal split component H with more than three bends. 
Then one of the two paths in the outer face of H has rotation less than —3 and we have the same 
situation as for a single edge, that is we can apply Lemma [3] to reduce the rotation of the opposite 
site and thus reduce the number of bends of H by one. Afterwards, we can set the flexibility of H 
down to the new number of bends ensuring that it is not increased later on. However, this only 
works if the negative rotation of the split component H lies in an inner face of G. On the outer 
face we can only increase to a rotation of —5 yielding an orthogonal representation such that every 
principal split component has at most three bends, or maybe four or five bends, if it has its negative 
rotation in the outer face. Note that this is essentially the same situation we also had in the proof 
of Theorem [2] In the following we show similarly that the number of bends can be reduced further, 
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Figure 9: The path between the new and the old root in the SPQR-tree containing \x (left). The 
whole graph G containing the principal split component H' corresponding to \x with respect to the 
new root and the principal split component H of the new root with respect to the old root (right). 

until either a unique innermost principal split component (where innermost means minimal with 
respect to inclusion) or the reference edge e re f may have more than three bends. 

First assume that e re f has more than three, that is four or five, bends and that there is a 
principal split component H with more than three bends having its negative rotation on the outer 
face. Let {s, t} be the corresponding split pair and let without loss of generality (t, s) be the path 
along H with rotation less than —3 where / is the outer face. Then the path 7Tf(s, t) contains the 
edge e re f = {u, v}, otherwise H would not be a principal split component. Moreover, rot(-Kf(t, s)) < 
—4 implies that rot(7Tj(s, t)) > —2 holds. As in the proof of Theorem [2] (compare with Figure[6^c)) 
the path nf(s,t) splits into the paths nf(s,u), nf(u,v) and Kf(v,t). Since itf(u,v) consists of the 
single edge e re f with more than three bends rot(7r/(it, v)) < —4 holds, implying that the rotation 
along ttj(s, u) or Ttf(v, t) is greater or equal to 0. This shows that 7iy(s, t) cannot be strictly directed 
from i to s and thus we can apply Lemma [3^1) to reduce the number of bends H has. Finally, 
there is no principal split component with more than three bends left and the reference edge 6 re f 
has at most five bends, which concludes this case. 

In the second case, e re f has at most three bends. We show that if there is more than one 
principal split component with more than three bends, then they hierarchically contain each other. 
Assume that the number of bends of no principal split component that has more than three bends 
can be reduced further. Assume further there are two principal split components Hi and H2 with 
respect to the split pairs {si, t±} and {S2, £2} that do not contain each other, that is without loss of 
generality the vertices t\, si, ti and S2 occur in this order around the outer face / when traversing 
it in counter-clockwise direction and 7T/(ii,si) and 7r/(i2, 82) belong to Hi and H2 respectively. 
Analogous to the case where e re f has more than three bends we can show that Lemma [3^1) can be 
applied to reduce the number of bends of Hi, which is a contradiction. Thus, either Hi is contained 
in H2 or the other way round. This shows that there is a unique principal split component H that 
is minimal with respect to inclusion having more than three bends. Due to the inclusion property, 
all nodes in the SPQR-tree corresponding to the principal split components with more than three 
bends lie on the path between the current root and the node corresponding to H. We denote the 
node corresponding to H by r and choose t to be the new root of the SPQR-tree T ■ Since the 
principal split components depend on the root chosen for T some split components may no longer 
be principal and some may become principal due to rerooting. Our claim is that all principal split 
components with more than three bends are no longer principal after rerooting and furthermore 
that all split components becoming principal can be enforced to have at most three bends. 

First note that the principal split component corresponding to a node [i in the SPQR-tree 
changes if and only if \x lies on the path between the old and the new root, that is between r 
and the Q-node corresponding to e re f. Since all principal split components (with respect to the 
old root) that have more than three bends also lie on this path, all these split components are 
no longer principal (with respect to the new root). It remains to deal with the new principal 
split components corresponding to the nodes on this path. Note that the new root r itself has no 
principal split component associated with it. Let / r be a node on the path between the new 
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and the old root and let H' be the new principal split component corresponding to \i with the 
poles s' and t! . Recall that H is the former principal split component corresponding to the new 
root r with the poles s and t. Note that H of course is still a split component, although it is not 
principal anymore. Figure [9] illustrates this situation. Now assume that H' has more than three 
bends. Then there are two possibilities, either it has its negative rotation on the outer face or in 
some inner face. If only the latter case arises we can easily reduce the number of bends down to 
three as we did before. In the remaining part of the proof we show that the former case cannot arise 
due to the assumption that the number of bends of H cannot be reduced anymore. Assume H' 
has its negative rotation in the outer face /, that is without loss of generality the path irf(t,s) 
belongs to H' and has rotation at most —4. Thus we have again the situation that the two split 
components H' and H both have a rotation of at most —4 in the outer face. Moreover, these two 
split components do not contain or overlap each other since s and t are not contained in H' as r 
is the new root and H does not contain s 1 or t' since fi is an ancestor of r with respect to the old 
root. Thus we could have reduced the number of bends of H before we changed the root, which is 
a contradiction to the assumption we made that the number of bends of principal split components 
with more than three bends cannot be reduced anymore. Hence, all new principal split components 
either have at most three bends or they have their negative rotation in some inner face. Finally, we 
obtain a valid orthogonal representation with at most three bends per principal split component 
with respect to r. □ 



5 Optimal Drawings with Fixed Planar Embedding 

All results from the previous sections deal with the case where we are only interested in the decision 
problem of whether a given graph has a valid drawing or not. More precisely, we always assumed to 
have a valid orthogonal representation of an instance of FlexDraw and showed that this implies 
that there exists another valid orthogonal representation with certain properties. In this section, 
we consider convex instances of the optimization problem OptimalFlexDraw. The following 
generic theorem shows that the results for FlexDraw that we presented so far can be extended 
to OptimalFlexDraw. 

Theorem 4. // the existence of a valid orthogonal representation of an instance of FlexDraw 
with positive flexibility implies the existence of a valid orthogonal representation with property P, 
then every convex instance of OptimalFlexDraw has an optimal drawing with property P. 

Proof. Let G be a convex instance of OptimalFlexDraw. Let further TZ be an optimal orthogonal 
representation. We can reinterpret G as an instance of FlexDraw with positive flexibility by 
setting the flexibility of an edge with p bends in TZ to max{p, 1}. Then TZ is obviously a valid 
orthogonal representation of G with respect to these flexibilities. Thus there exists another valid 
orthogonal representation TZ' having property P. It remains to show that cost (72.') < cost (72.) holds 
when going back to the optimization problem OptimalFlexDraw. However, this is clear for the 
following reason. Every edge e has as most as many bends in TZ' as in TZ except for the case where e 
has one bend in TZ' and zero bends in TZ. In the former case the monotony of cost e (-) implies that 
the cost did not increase. In the latter case e causes the same amount of cost in TZ as in TZ' since 
cost e (0) = cost e (l) = b e holds for convex instances of OptimalFlexDraw. Note that this proof 
still works, if the cost functions are only monotone but not convex. □ 

It follows that every convex 4-planar graph has an optimal drawing that is nice since Theorem [4] 
shows that Theorem [3] can be applied. Thus, it is sufficient to consider only nice drawings when 
searching for an optimal solution, as there exists a nice optimal solution. This is a fact that we 
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Figure 10: Split components with as few bends as possible. 



crucially exploit in the next section since although the cost function of a principal split component 
may be non-convex, we can show that it is convex in the interval that is of interest when only 
considering nice drawings. 



6 Optimal Drawings with Variable Planar Embedding 

All results we presented so far were based on a fixed planar embedding of the input graph G. In 
this section we present an algorithm that computes an optimal drawing of G in polynomial time, 
optimizing over all planar embeddings of G. Our algorithm crucially relies on the existence of a 



nice drawing among all optimal drawings of G. For biconnected graphs (Section 6.1) we present a 
dynamic program that computes the cost function of all principal split components bottom-up in 
the SPQR-tree with respect to a chosen root. To compute the optimal drawing among all drawings 
that are nice with respect to the chosen root, it remains to consider the embeddings of the root itself. 
If we choose every node to be the root once, this directly yields an optimal drawing of G taking 



all planar embeddings into account. In Section 6.2 we extend our results to connected graphs that 



are not necessarily biconnected. To this end we first modify the algorithm for biconnected graphs 
such that it can compute an optimal drawing with the additional requirement that a specific vertex 
lies on the outer face. Then we can use the BC-tree to solve OptimalFlexDraw for connected 
graphs. We use the computation of a minimum-cost flow in a network of size re as a subroutine and 



denote the consumed running time by Tfl ow (n). In Section 6.3 we consider which running time we 
actually need. 



6.1 Biconnected Graphs 

In this section we always assume G to be a biconnected 4-planar graph forming a convex instance 
of OptimalFlexDraw. Let T be the SPQR-tree of G. As defined before, an orthogonal represen- 
tation is optimal if it has the smallest possible cost. We call an orthogonal representation r-optimal 
if it has the smallest possible cost among all orthogonal representation that are nice with respect to 
the root r. We say that it is (t,£) -optimal if it causes the smallest possible amount of cost among 
all orthogonal representations that are nice with respect to r and induce the planar embedding £ on 
skel(7~). In this section we concentrate on finding a (t, £ )-optrmal orthogonal representation with 
respect to a root r and a given planar embedding £ of skel(r). Then a r-optimal representation 
can be computed by choosing every possible embedding of skel(r). An optimal solution can then 
be computed by choosing every node in T to be the root once. 

In Section [4] we extended the terms "bends" and "flexibility", which were originally defined for 
single edges, to arbitrary principal split components with respect to the chosen root. We start out by 
making precise what we mean with the cost function cost#(-) of a principal split component H with 
poles s and t. Recall that the number of bends of H with respect to an orthogonal representation 
S with s and t on the outer face / is defined to be max{| rots(7r/(s, t))\,\ rots(7T/-(i, s))\}. Assume 
S is the nice orthogonal representation of H that has the smallest possible cost among all nice 
orthogonal representations with p bends. Then we essentially define cost h(p) to be the cost of 



19 



Figure 11: A single vertex can be replaced by a split component with three bends. 



S. However, with this definition the cost function of H is not defined for all p G No since H 
does not have an orthogonal representation with zero bends at all, if deg(s) > 1 or deg(i) > 1, as 
at least one of the paths nf(s,t) and nf(t,s) has negative rotation in this case. More precisely, 
if deg(s) + deg(i) > 2, then H has at least one bend, and if deg(s) + deg(t) > 4, then H has 



at least two bends. Figure 10 shows for each combination of degrees a small example with the 
smallest possible number of bends. In these two cases we formally set costff (0) = costff (1) and 
cost # (0), cost if (1) = cost h (2), respectively. Thus, we only need to compute the cost functions for at 
least [(deg(s)+deg(i)— 2)/2] bends. We denote this lower bound by In = [(deg(s) +deg(t) — 2)/2~|. 
Hence, it remains to compute the cost function costff (p) for p S [£h, 3]. For more than three bends 
we formally set the cost to oo. Note that the definition of the cost function only considers nice 
orthogonal representations (including that they are tight). As a result of this restriction the cost for 
an orthogonal representation with p bends might be less than costff (p). However, due to Theorem [3] 
in combination with Theorem [4] we know that optimizing over nice orthogonal representations is 
sufficient to find an optimal solution. 

As for single edges, we define the base cost bn of the principal split component H to be costff (0). 
We will see that the cost function costff (•) is monotone and even convex in the interval [0, 3] (except 
for a special case) and thus the base cost is the smallest possible amount of cost that has to be 
payed for every orthogonal drawing of H. The only exception is the case where deg(s) = deg(i) = 3. 
In this case H has at least two bends and thus the cost function costff (•) needs to be considered 
only on the interval [2,3]. However, it may happen that costff (2) > costff (3) holds in this case. 
Then we set the base cost bn to costff (3) such that the base cost bn is really the smallest possible 
amount of cost that need to be payed for every orthogonal representation of H. We obtain the 
following theorem. 

Theorem 5. // the poles of a principal split component do not both have degree 3, then its cost 
function is convex on the interval [0,3]. 

Before showing Theorem [5] we just assume that it holds and moreover we assume that the cost 
function of every principal split component is already computed. We first show how these cost 
functions can then be used to compute an optimal drawing. To this end, we define a flow network 
on the skeleton of the root r of the SPQR-tree, similar to Tamassias flow network [15 . The cost 
functions computed for the children of r will be used as cost functions on arcs in the flow network. 
As we can only solve flow networks with convex costs we somehow have to deal with potentially 
non-convex cost functions for the case that both endvertices of a virtual edge have degree 3 in its 
expansion graph. Our strategy is to simply ignore these subgraphs by contracting them into single 
vertices. Note that the resulting vertices have degree 2 since the poles of graphs with non-convex 
cost functions have degree 3. The process of replacing the single vertex in the resulting drawing by 
the contracted component is illustrated in Figure [TTJ The following lemma justifies this strategy. 



Lemma 4. Let G be a biconnected convex instance of OptimalFlexDraw with t -optimal orthog- 
onal representation 7Z and let H be a principal split component with non-convex cost function and 
base cost bn- Let further G' be the graph obtained from G by contracting H into a single vertex 
and let 1Z' be a r-optimal orthogonal representation of G' . Then cost(7£) = cost(JZ') + bn holds. 
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Figure 12: (a) The structure of the flow network N for the case that r is an R-node with skel(7~) = 
K4. The outer face is split into several gray boxes to improve readability, (b) A flow together with 
the corresponding orthogonal representation. The numbers indicate the amount of flow on the arcs. 
Undirected edges imply flow, directed arcs without a number have flow 1. 



Proof. Assume we have a r-optimal orthogonal representation 1Z of G inducing the orthogonal 
representation S on H. As H has either two or three bends we can simply contract it yielding 
an orthogonal representation 1Z' of G with cost(7£') = cost(7£) — cost(5) < cost(7£) — bu- The 
opposite direction is more complicated. Assume we have an orthogonal representation 1Z' of G' , 
then we want to construct an orthogonal representation 1Z of G with cost(7£) = cost(7£/) + bu- 
Let S be an orthogonal representation of H causing only bu cost. Since cost#(-) was assumed to 
be non-convex, S needs to have three bends. It is easy to see that 1Z' and S (or S' obtained from S 
by mirroring the drawing) can be combined to an orthogonal representation of G if the two edges 
incident to the vertex v in G' corresponding to H have an angle of 90° between them. However, this 
can always be ensured without increasing the costs of 1Z' . Let e\ and e 2 be the edges incident to v 
and assume they have an angle of 180° between them in both faces incident to v. If neither e\ nor e 2 
has a bend, the flex graph contains the cycle around v due to the fact that e\ and e 2 have positive 
flexibilities. Bending along this cycles introduces a bend to each of the edges, thus we can assume 
without loss of generality that e\ has a bend in 1Z' . Moving v along the edge e\ until it reaches this 
bend decreases the number of bends on e\ by one and ensures that v has an angle of 90° in one of its 
incident faces. Thus we can replace v by the split component H with orthogonal representation S 
having cost bn yielding an orthogonal representation 1Z of G with cost (TZ) = cost(7£') + bn- □ 

When computing a (r, £)-optimal orthogonal representation of G we make use of Lemma [4] in 
the following way. If the expansion graph H corresponding to a virtual edge e in skel(r) has a non- 
convex cost function, we simply contract this virtual edge in skel(r). Note that this is equivalent 
to contracting H in G. We can then make use of the fact that all remaining expansion graphs have 
convex cost functions to compute a (r, £ )-optimal orthogonal representation of the resulting graph 
yielding a (r, £)-optimal orthogonal representation of the original graph G since the contracted 
expansion graphs can be inserted due to Lemma |4| Note that expansion graphs with non convex 
cost functions can only appear if the root is a Q- or an S-node. In the skeletons of P- and R-nodes 
every vertex has degree at least three, thus the poles of an expansion graph cannot have degree 3 
since G has maximum degree 4. 

Now we are ready to define the flow network iV^ on skel(r) with respect to the fixed embedding £ 
of skel(r); see Figure [l2fa) for an example. For each vertex v, each virtual edge e and each face / 
in skel(r) the flow network iV^ contains the nodes v, e and /, called vertex node, edge node and 
face node, respectively. The network contains the arcs (v,f) and (/, v) with capacity 1, called 
vertex-face arcs, if the vertex v and the face / are incident in skel(r). For every virtual edge e we 
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add edge-face arcs (e, /) and (/, e), if / is incident to e. We use cost#(-) — as cost function 
of the arc (f,s), where H is the expansion graph of the virtual edge e. The edge- face arcs (e, /) 
in the opposite direction have infinite capacity with cost. It remains to define the demand of 
every node in N £ . Every inner face has a demand of 4, the outer face has a demand of —4. 
An edge node e stemming from the edge e = {s, t} with expansion graph H has a demand of 
deg^(s) + deg^(i) — 2, where deg H (v) denotes the degree of v in H. The demand of a vertex 
node v is 4 - deg G (t;) - deg skel(r) (v). 

In the flow network N £ the flow entering a face node / using a vertex-face arc or an edge- 
face arc is interpreted as the rotation at the corresponding vertex or along the path between the 



poles of the corresponding child, respectively; see Figure 12 ^b) for an example. Incoming flow is 
positive rotation and outgoing flow negative rotation. Let bn 1 bu k be the base costs of the 
expansion graphs corresponding to virtual edges in skel(r). We define the total base costs of r to 
be b T = J2i °Hi ■ Note that the total base costs of r are a lower bound for the costs that have to be 
paid for every orthogonal representation of G. We show that an optimal flow in N £ corresponds 
to a (r, £)-optimal orthogonal representation 1Z of G. Since the base costs do not appear in the 
flow network, the costs of the flow and its corresponding orthogonal representation differ by the 
total base costs b T , that is cost (JZ) = cost(0) + b T . We obtain the following lemma. 

Lemma 5. Let G be a biconnected convex instance of OptimalFlexDraw, let T be its SPQR- 
tree with root t and let £ be an embedding of skel(r) . If the cost function of every principal split 
component is known, a (r,£)-optimal solution can be computed in 0(Tfl ow (| skel(r)|)) time. 

Proof. As mentioned before, we want to use the flow network N £ to compute an optimal orthogonal 
representation. To this end we show two directions. First, given a (r, £)-optimal orthogonal 
representation 7Z, we obtain a feasible flow <f> in N £ such that cost(^) = cost {JZ) — b T , where b T 
are the total base costs. Conversely, given an optimal flow <fi in N , we show how to construct an 
orthogonal representation 1Z such that cost (JZ) = cost(0) + b T . As the flow network N £ has size 
0(| skel(r)|), the claimed running time follows immediately. 

Let 1Z be a (r, £ )-optimal orthogonal representation of G. As we only consider nice and thus 
only tight drawings we can assume the orthogonal representation 1Z to be tight. Recall that being 
tight implies that the poles of the expansion graph of every virtual edge have a rotation of 1 in 
the internal faces. We first show how to assign flow to the arcs in N £ . It can then be shown that 
the resulting flow is feasible and causes cost(7£) — b T cost. For every pair of vertex-face arcs (/, v ) 
and (v, f) in N £ there exists a corresponding face / in the orthogonal representation IZolG and we 
set 4>((v , /)) = rot(vf). Let e = {s, t} be a virtual edge in skel(/x) incident to the two faces fi and fi- 
Without loss of generality let nf 1 (s,t) be the path belonging to the expansion graph of e. Then 
irf 2 (t,s) also belongs to H. We set 4>((e,fi)) = rot7^.(7T/ 1 (s, t)) and /2)) = r ot7e(7T/ 2 (t, s)). For 
the resulting flow (j) we need to show that the capacity of every arc is respected, that the demand 
of every vertex is satisfied, and that cost(0) = cost(7£) — b T holds. 

First note that the flow on the vertex-face arcs does not exceed the capacities of 1 since every 
vertex has degree at least 2. Since no other arc has a capacity, it remains to deal with the demands 
and the costs. 

For the demands we consider each vertex type separately. Let / be a face node. The total 
incoming flow entering / is obviously equal to the rotation in 1Z around the face /. As 1Z is an 
orthogonal representation this rotation equals to 4 (—4 for the outer face), which is exactly the 
demand of /. Let e be an edge node corresponding to the expansion graph H with poles s and t. 
Recall that dem(e) = deg^(s) + deg^(i) — 2 is the demand of e. Figure [i~3|a) illustrates the 
demand of a virtual edge. Let S be the orthogonal representation induced on H by 1Z and let / 
be the outer face of S. Clearly, the flow leaving e is equal to rot7j,(7r/ 1 (s, t)) + rot-^(7rj 2 (t, s)) = 
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Figure 13: (a) Illustration of the demand of virtual edges, (b) Rotation of poles in the outer face, 
depending on the degree. 



rot s (tt f(s,t)) + rots(irf(t, s)). Since / is the outer face of H, the total rotation around this faces 
sums up to —4. The rotation of the pole s in the outer face / is deg^(s) — 3, see Figures [l3^b) , and 
the same holds for t. Thus we have rots{^f{s, t)) +rot l s(7r s)) +deg^(s) — 3 + deg^(i) — 3 = —4. 
This yields for the outgoing flow rot,s(7r/(s, t)) + rots(7r/(t, s)) = 2 — deg H (s) — deg H (t), which is 
exactly the negative demand of e. It remains to consider the vertex nodes. Let v be a vertex node, 
recall that dem(u) = 4 — deg G (u) — deg ske i( T )(t>) holds. The outgoing flow leaving v is equal to 
the summed rotation of v in faces not belonging to expansion graphs of virtual edges in skel(r). 
As 1Z is an orthogonal representation, the total rotation around every vertex v is 2 • (deg G (v) — 2). 
Moreover, v is incident to deg sko i( r )(v) faces that are not contained in expansion graphs of virtual 
edges of skel("r). Thus there are deg G (?;) — deg skcl ( T )(v) faces incident to v belonging to expansion 
graphs. As we assumed that the orthogonal representation of every expansion graph is tight, the 
rotation of v in each of these faces is 1 . Thus the rotation of v in the remaining faces not belonging 
to expansion graphs is 2 • (deg G (v) — 2) — (deg G (v) — deg skel ( T ) {y )) . Rearrangement yields a rotation, 
and thus an outgoing flow, of deg G {v) + deg ske i( T )(-u) — 4, which is the negative demand of v. 

To show that cost(^) = cost (7?.) — b T holds it suffices to consider the flow on the edge- face arcs as 
no other arcs cause cost. Let e be a virtual edge and let /i and ji the two incident faces. The flow 
entering f\ or /2 does not cause any cost, as (e, /i) and (e, $2) have infinite capacity with cost. 
Thus only flow entering e over the arcs (fi,e) and (/2,e) may cause cost. Assume without loss of 
generality that the number of bends p the expansion graph H of e has is determined by the rotation 
along nf^Sjt), that is p = — votn^^ (s, t)). Let p' = — mt-ji^ f 2 {t, s)) be the negative rotation 
along the path 7Tj 2 (i, s) in the face /2- Note that </>((fi,s)) = p and (/>((/2,e)) = p' . Obviously, the 
flow on (fx,s) causes the cost cost#(p) — bu- We show that the cost caused by the flow on (/2,e) 
is 0. If p' < this is obviously true, as there is no flow on the edge (/2,^). Otherwise, < p' < p 
holds. It follows that the smallest possible number of bends In every orthogonal representation 
of H has lies between p' and p. It follows from the definition of cost#(-) and from the fact that 
all cost functions are convex that cost h(p') = bn- To sum up, the total cost on edge- face arcs 
incident to the virtual edge e is equal to the cost caused by its expansion graph H with respect to 
the orthogonal representation 1Z minus the base cost bn- As neither (ft nor 1Z have additional cost 
we obtain cost((/>) = cost (1Z) — b T . 

It remains to show the opposite direction, that is given an optimal flow <j> in iV^, we can 
construct an orthogonal representation 1Z of G such that cost (JZ) = cost(i^) + b T . This can be done 
by reversing the construction above. The flow on edge-face arcs determines the number of bends for 
the expansion graphs of each virtual edge. The cost functions of these expansion graphs guarantee 
the existence of orthogonal representations with the desired rotations along the paths between the 
poles, thus we can assume to have orthogonal representations for all children. We combine these 
orthogonal representations by setting the rotations between them at common poles as specified by 
the flow on vertex- face arcs. It can be easily verified that this yields an orthogonal representation 
of the whole graph G by applying the above computation in the opposite direction. □ 

The above results rely on the fact that the cost functions of principal split components are 
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convex as stated in Theorem [5] and that they can be computed efficiently. In the following we show 
that Theorem [5] really holds with the help of a structural induction over the SPQR-tree. More 
precisely, the cost functions of principal split components corresponding to the leaves of T are the 
cost functions of the edges and thus they are convex. For an inner node p we assume that the 
pertinent graphs of the children of p have convex cost functions and show that H = pert(^) itself 
also has a convex cost function. The proof is constructive in the sense that it directly yields an 
algorithm to compute these cost functions bottom up in the SPQR-tree. 

Note that we can again apply Lemma [4] in the case that the cost function of the expansion 
graph of one of the virtual edges in skel(yu) is not convex due to the fact that both of its poles 
have degree 3. This means that we can simply contract such a virtual edge (corresponding to 
a contraction of the expansion graph in H), compute the cost function for the remaining graph 
instead of H and plug the contracted expansion graph into the resulting orthogonal representations. 
Thus we can assume that the cost function of each of the expansion graphs is convex, without any 
exceptions. 

The flow network N £ that was introduced to compute an optimal orthogonal representation 
in the root of the SPQR-tree can be adapted to compute the cost function of the principal split 
component H corresponding to a non-root node p. To this end we have to deal with the parent 
edge, which does not occur in the root of T, and we consider a parameterization of A^ to compute 
several optimal orthogonal representations with a prescribed number of bends, depending on the 
parameter in the flow network. Before we describe the changes in the flow network we need to make 
some considerations about the cost function. By the definition of the cost function it explicitely 
optimizes over all planar embeddings of skel(yu). Moreover, as the cost function cost h{p) depends 
on the number of bends p a graph H has, it implicitly allows to flip the embedding of H since 
the number of bends is defined as max{| rot(7r(s, t))\, | rot(7r(i, s))|}. However, the flow network TV 
can only be used to compute the cost function for a fixed embedding. Thus we define the partial 
cost function costfj (p) of H with respect to the planar embedding £ of skel(p) to be the smallest 
possible cost of an orthogonal representation inducing the planar embedding £ on skel(/i) with p 
bends such that the number of bends is determined by -Kf(s,t), that is rot(7T/(s, t)) = —p, where / 
is the outer face. Note that the minimum over the partial cost functions cost^(-) and cost^(-), 
where £' is obtained by flipping the embedding £ of skel(/i) yields a function describing the costs 
of H with respect to the embedding £ of skel(/i) depending on the number of bends H has (and 
not on the rotation along vrj(s, t) as the partial cost function does). Obviously, minimizing over all 
partial cost functions yields the cost function of H . 

The flow network N £ is defined as before with the following modifications. The parent edge of 
skel(^t) does not have a corresponding edge node. Let f\ and /2 be the faces in skel(/x) incident 
to the parent edge. These two faces together form the outer face / of H, thus we could merge 
them into a single face node. However, not merging them has the advantage that the incoming 
flow in fi and f2 corresponds to the rotations along irf(s,t) and 7T/(i, s), respectively (it might be 
the other way round but we can assume this situation without loss of generality). Thus, we do not 
merge f\ and /2, which enables us to control the number of bends of H by setting the demands 
of fi and /2- This is also the reason why we remove the vertex- face arcs between the poles and the 
two faces f\ and /2- Before we describe how to set the demands of fi and /2, we fit the demands 
of the poles to the new situation. As we only consider tight orthogonal representations we know 
that the rotation at the poles s and t in all inner faces is 1. Thus, we set dem(s) = 2 — deg skel (^)(s) 
and dem(i) = 2 — deg skel ( M ) (t) as this is the number of faces incident to s and t, respectively, after 
removing the vertex- face arcs to /i and /2- With these modifications the only flow entering f\ and /2 
comes from the paths vrj(s, t) and 71/ (t, s), respectively. As the total rotation around the outer face 
is —4 and the rotation at the vertices s and t is deg H (s) — 3 and deg^(i) — 3, respectively, we have to 
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ensure that dem(/i) + dem(/2) = 2 — deg^(s) — deg H (t). As mentioned before, we assume without 
loss of generality that ixt (s, t) belongs to the face /i and 7r/ (i, s) belongs to f2- Then the incoming 
flow entering f% corresponds to rot(-7r/(s, t)) of an orthogonal representation. We parameterize N £ 
with respect to the faces fi and fi starting with dem(/i) = and dem(/2) = 2 — deg^(s) — deg H (t). 
It obviously follows that an optimal flow in N £ with respect to the parameter p corresponds to 
an optimal orthogonal representation of H that induces £ on skel(^) and has a rotation of —p 
along 7Tf(s,t). Thus, up to the total base costs b^, the cost function of the flow network equals 
to the partial cost function of H on the interval [-£#,3], that is cost N e(p) + 6 M = cost^(p) for 
£h < p < 3. To obtain the following lemma it remains to show two things for the case that 
deg(s) + deg(i) < 6. First, cost A rf(p) and thus each partial cost function is convex for £n < P < 3. 
Second, the minimum over these partial cost functions is convex. 

Lemma 6. // Theorem^ holds for each principal split component corresponding to a child of the 
node p in the SPQR-tree, then it also holds for pert(/i). 

Proof. As mentioned before, we can use the flow network N £ to compute the partial cost function 
costff(p) for tu < p < 3 since cost^(p) = cost N e(p) + b^ holds on this interval. In the following 
we only consider the case where degj^(s) + deg^(i) < 6 holds for the poles s and t. For the case 
deg^(s) = degn(t) = 3 we do not need to show anything. To show that the partial cost function is 
convex we do the following. First, we show that cost^(p) is minimal for p = £n- This implies that 
the cost function cost N e(p) of the flow network is minimal for p = po < £n- Then Theorem [T] can 
be applied showing that cost N e(p) is convex for p E [po,oo] yielding that the partial cost function 
costf^/}) is convex for p E [£h, 3]. Thus, it remains to show that cost^(p) is minimal for p = £n to 
obtain convexity for the partial cost functions. 

Let S be an orthogonal representation of H with p E [£n, 3] bends such that 7r/(s, t) determines 
the number of bends, that is rots (71V (s,t)) = — Pi where / is the outer face of H. We show the 
existence of an orthogonal representation S' with rots/(7r/(s, i)) = —in and cost (5') < cost(5). 
Since we assume S to be tight, the rotations at the poles rots(sj) and rots(tf) only depend on the 
degree of s and t. More precisely, we have rots(s/) = deg^(s) — 3 and the same holds for t. Since 
the total rotation around the outer face / is —4 the following equation holds. 

rot lS (7r / (t, s)) = p + 2- deg H (s) - deg H (t) (1) 

In the following we show that Tots(^f(t, s)) > holds if the number of bends p exceeds £jj- Then 
Corollary [T] in combination with Theorem [4] can be used to reduce the rotation along nf(t,s) 
and thus reduce the number of bends by 1, yielding finally an orthogonal representation with in 
bends determined by irf(s,t). Recall that the lower bound for the number of bends was defined 
as In = [(deg(s) + deg(t) — 2)/2~|. First consider the case that deg^(s) + deg^(i) is even (and of 
course less than 6). Then Equation ([T]) yields rots(7r/(f, s)) = p — 2£n- If p is greater than £n this 
yields iots(^f(t, s)) > —£n- Since £n is at most 1 in the case that deg(s) + deg(i) is even and less 
than 6, this yields rot,s(7r/(t, s)) > —1. The case that deg^(s) + deg H (t) is odd works similarly. 
Then Equation ([TJ yields rots(7r/(£, s)) = p — 2£jj + l. As before p is assumed to be greater than £n 
yielding vots(TTf(t, s)) > —in + 1. As £n is at most 2 we again obtain rots(7T/(t, s)) > —1, which 
concludes the proof that the partial cost functions are convex. 

It remains to show that the minimum over the partial cost functions is convex. First assume 
that p is an R-node. Then its skeleton has only two embeddings £ and £' where £' is obtained by 
flipping £. We have to show that the minimum over the two partial cost functions cost^(-) and 
costff(-) remains convex. For the case that deg(s) +deg(t) = 5 the equation £n = 2 holds and thus 
we only have to show convexity on the interval [2, 3]. Obviously, cost^(-) is convex on this interval 



25 



if and only if cost#(2) < cost#(3). As this is the case for both partial cost functions, it is also 
true for the minimum. For deg(s) + deg(i) < 5 we first show that cost#(^#) = costf^(£#) holds. 
For the case that deg(s) + deg(t) is even this is clear since mirroring an orthogonal representa- 
tion S with rots(7T/(s, t)) = —in inducing £ on skel(//) yields an orthogonal representation S' with 
rots'(7r/(s, t)) = —i-H inducing £' on skel(/i). For the case that deg(s) +deg(i) = 3, the orthogonal 
representation S with rotation —1 along irf(s,t) can also be mirrored yielding S' with rotation 
along iTf(s,t). By Corollary [T] this rotation can be reduced to —1 without causing any additional 
cost. As this construction also works in the opposite direction we have cost^(fe) = cost jj(£h) 
for all cases. Moreover, cost^(O) = cost^(l) holds by definition, if deg(s) + deg(i) > 2. If 
deg(s) = deg(t) = 1 this equation is also true as the rotation along irf(s,t) of an orthogonal 
representation can be reduced by 1 if it is 0, again due to Corollary [TJ Thus it remains to show 
that the cost function cost#(-) defined as the minimum of cost^(-) and cost^(-) is convex on the 
interval [1,3]. 

Assume for a contradiction that cost h(p) is not convex for p £ [1,3], that is A costal) > 
Acost#(2). Assume without loss of generality that cost#(3) = cost^(3) holds. As we showed before 
cost#(l) = costf^(l) also holds. Since cost#(2) is the minimum over cost^(2) and cost^(2) we 
additionally have cost#(2) < costf^(2). This implies that the inequalities Acost^(l) > Acost#(l) 
and Acost#(2) < Acost#(2) hold, yielding that the partial cost function cost^(p) is not convex 
for p £ [1,3], which is a contradiction. Thus cost#(-) is convex. 

The case that p is a P-node works similar to the case that p is an R-node. If p has only two 
children, its skeleton has only two embeddings £ and £' obtained from one another by flipping. 
Thus the same argument as for R-nodes applies. If p has three children, then deg(s) = deg(t) = 3 
holds and thus we do not have to show convexity. Note that in the case deg(s) = deg(t) = 3 the 
resulting cost function can be computed by taking the minimum over the partial cost functions with 
respect to all embeddings of skel(/x), although it may by non-convex. If p is an S-node, we have 
a unique embedding and thus the partial cost function with respect to this embedding is already 
the cost function of H. Note that considering only the rotation along nf(s,t) for the partial cost 
function is not a restriction, as S-nodes are completely symmetric. □ 

Lemma [6] together with the fact that the cost function of every edge is convex shows that 
Theorem [5] holds, that is the cost functions of all principal split components are convex on the 
interesting interval [0, 3] except for the special case where both poles have degree 3. However, 
this special case is easy to handle as principal split components of this type with non-convex 
cost functions can be simply contracted to a single vertex by Lemma [4] Moreover, the proof is 
constructive in the sense that it shows how the cost functions can be computed efficiently bottom 
up in the SPQR-tree. For each node p we have to solve a constant number of minimum-cost flow 
problems in a flow network of size 0{\ skel(/i)|). As the total size of all skeletons in T is linear 
in the number n of vertices in G, we obtain an overall C(Tfl ow (n)) running time to compute the 
cost functions with respect to the root r. Finally, Lemma [5] can be applied to compute an optimal 
orthogonal representation with respect to a fixed root and a fixed embedding of the root's skeleton 
in 0(Tfl ow (| skel(r)|)) time. To compute an overall optimal solution, we have to compute a (r, £ )- 
optimal solution for every root r and every embedding £ of skel(r). The number of embeddings of 
skel(r) is linear in the size of skel(r) (since P-nodes have at most degree 4) and the total size of all 
skeletons is linear in n. We obtain the following theorem. 

Theorem 6. OptimalFlexDraw can be solved in 0(n ■ Jfl ow (n)) time for convex biconnected 
instances. 
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6.2 Connected Graphs 



In this section we extend the result obtained in Section 6.1 to the case that the input graph G 
contains cutvertices. Let B be the BC-tree of G rooted at some B-node j3. Then every Block except 
for (5 has a unique cutvertex as parent and we need to find optimal orthogonal representations with 
the restriction that this cutvertex lies on the outer face. We claim that we can then combine these 
orthogonal representations of the blocks without additional cost. 

Unfortunately, with the so far presented results we cannot compute the optimal orthogonal 
representation of a biconnected graph considering only embeddings where a specific vertex v lies 
on the outer face. We may restrict the embeddings of the skeletons we consider when traversing 
the SPQR-tree bottom up to those who have v on the outer face. However, we can then no longer 
assume that the cost functions we obtain are symmetric. To deal with this problem, we present a 
modification of the SPQR-tree, that can be used to represent exactly the planar embeddings that 
have v on the outer face and are represented by the SPQR-tree rooted at a node r. 

Let r be the root of the SPQR-tree T ■ If v is a vertex of skel(r) , then restricting the embeddings 
of skel(r) to those who have v on the outer face of skel(r) forces v to be on the outer face of the 
resulting embedding of G. Otherwise, v is contained in the expansion graph of a unique virtual 
edge s in skel(r), we say that v is contained in e. Obviously, e has to be on the outer face of the 
embedding of skel(r). However, this is not sufficient and it depends on the child fi of r corresponding 
to e whether v lies on the outer face of the resulting embedding of G. Let £ T be an embedding 
of skel(r) having e on the outer face and let s and t be the endpoints of e. Then there are two 
possibilities, either e = {s,t} has the outer face to the left or to the right, where the terms "left" 
and "right" are with respect to an orientation from t to s. Assume without loss of generality that 
the outer face lies to the right of e and consider the child \x of r corresponding to e. As T is 
rooted, we consider only embeddings of skel(/u) that have the parent edge {s,t} on the outer face. 
As the choice of the outer face of skel(fi) does not have any effect on the resulting embedding, we 
can assume that {s,t} lies to the left of skel(^), that is the inner face incident to {s,t} lies to the 
right of {s, t} with respect to an orientation from t to s. A vertex contained in skel(/i) then lies 
obviously on the outer face of the resulting embedding of G if and only if it lies on the outer face 
of the embedding of skel(fi). Thus, if v is contained in skel^), restricting the embedding choices 
such that v lies on the outer face of skel(/x) forces v to be on the outer face of G. Note that in this 
case \i is either an R- or an S-node. For S-nodes there is no embedding choice and every vertex in 
skel(^t) lies on the outer face in this embedding. If [i is an R-node, there are only two embeddings 
and either v lies on the outer face of exactly one of them or in none of them. In the latter case the 
SPQR-tree with respect to the root r does not represent an embedding of G with v on the outer 
face at all. 

Assume that v is not contained in skel(/x). Then it is again contained in a single virtual edge e' 
and it is necessary that e' lies on the outer face of the embedding of skel(/u). Moreover, it depends 
on the child of [i corresponding to e' whether v really lies on the outer face. Note that fixing e' 
on the outer face completely determines the embedding of skel(/i) if it is not a P-node. If fj, is a 
P-node, the virtual edge e' has to be the rightmost, whereas the order of all other virtual edges can 
be chosen arbitrarily. If this is the case we split the P-node into two parts, one representing the 
fixed embedding of e', the other representing the choices for the remaining edges; see Figure llffia) . 
More precisely, we split \i into two P-nodes, the first one containing the parent edge {s,t}, the 
edge e' and a new virtual edge corresponding to the second P-node, which is inserted as child. The 
skeleton of the second P-node contains a parent edge corresponding to the first P-node and the 
remaining virtual edges that were contained in skel(^) but are not contained in the first P-node. 
The children of \i are attached to the two P-nodes depending on where the corresponding virtual 
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Figure 14: (a) Splitting a P-node into two P-nodes, the vertex v fixed to the outer face is contained 
in the thick edges, (b) Contracting the path from the root to the node containing v in its skeleton. 



edges are. Note that by splitting the P-node //, the virtual edge e 1 can no longer be in between two 
other virtual edges in fi. However, this is a required restriction, thus we do not loose embeddings 
that we want to represent. Moreover, the new P-node containing the virtual edge e' that need 
to be fixed to the outer face contains only two virtual edges (plus the parent edge) and thus the 
embedding of its skeleton is completely fixed by requiring e' to be on the outer face. 

To sum up, if skel(r) contains v, then we simply have to choose an embedding of skel(r) with v 
on the outer face. Otherwise, we have to fix the virtual edge containing v to the outer face and 
additionally have to consider the child of r corresponding to this virtual edge. For the child we 
then have essentially the same situation. Either v is contained in its skeleton, then the embedding 
is fixed to the unique embedding having v on the outer face or v is contained in some virtual 
edge. However, then the embedding of the skeleton is again completely fixed (P-nodes have to be 
split up first) and we can continue with the child corresponding to the virtual edge containing v. 
This yields a path of nodes starting with the root r having a completely fixed embedding only 
depending on the embedding £ T chosen for skel(r). As the nodes on the path do not represent any 
embedding choices, we can simply contract the whole path into a single new root node, merging 
the skeletons on the path, such that the embedding of the new skeleton of the root is still fixed. 
This contraction is illustrated in Figure 14 b). More precisely, let r be the root and let e be the 



edge containing v, corresponding to the child \i. Then we merge r and /i by replacing e in r by the 
skeleton of m without the parent edge. The children of m are of course attached to the new root 
t' since skel(r') contains the corresponding virtual edges. As mentioned before, the embedding of 
skel(/i) was fixed by the requirement that v is on the outer face, thus the new skeleton skel(r') has 
a unique embedding £ T i inducing £ T on skel(r) and having v or the new virtual edge containing v 
on the outer face. The procedure of merging the root with the child corresponding to the virtual 
edge containing v is repeated until v is contained in the skeleton of the root. We call the resulting 
tree the restricted SPQR-tree with respect to the vertex v and to the embedding £ T of the root. 

To come back to the problem OptimalFlexDraw, we can easily apply the algorithm presented 
in Section 6.1 to the restricted SPQR-tree. All nodes apart from the root are still S-, P-, Q- or 



R-nodes and thus the cost functions with respect to the corresponding pertinent graphs can be 
computed bottom up. The root r may have a more complicated skeleton, however, its embedding 
is fixed, thus we can apply the flow algorithm as before, yielding an optimal drawing with respect 
to the chosen root r and to the embedding £ T of skel(r) with the additional requirement that v 
lies on the outer face. Since the restricted SPQR-tree can be easily computed in linear time for 
a chosen root r and a fixed embedding £ of skel(r), we can compute a (r, £ )-optimal orthogonal 
representation with the additional requirement that v lies on the outer face in Tfj ow (n) time, yielding 
the following theorem. 
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Theorem 7. OptimalFlexDraw with the additional requirement that a specific vertex lies on 
the outer face can be solved in 0{n ■ Tfl ow (n)) time for convex biconnected instances. 

As motivated before, we can use the BC-tree to solve OptimalFlexDraw for instances that 
are not necessarily biconnected. We obtain the following theorem. 

Theorem 8. OptimalFlexDraw can be solved in 0(n 2 • Tfl ow (n)) time for convex instances. 

Proof. Let G be a convex instance with positive flexibility of OptimalFlexDraw and let B be its 
BC-tree rooted at some B-node f5. We show how to find an optimal drawing of G, optimizing over all 
embeddings represented by B with respect to the root (3. Then we can simply choose every B-node 
in B to be the root once, solving OptimalFlexDraw. The algorithm consumes 0(n-Tfl ow (n)) time 
for each root (3 and thus the overall running time is 0(n 2 ■ Ja ow (n)). For the block corresponding 
to the root (3 we use Theorem [6] to find the optimal orthogonal representation. For all other blocks 
we use Theorem [7] to find the optimal orthogonal representation with the cutvertex corresponding 
to the parent in B on the outer face. It remains to stack these orthogonal representations together 
without causing additional cost. This can be easily done, if a cutvertex that is forced to lie on the 
outer face has all free incidences in the outer face and every other cutvertex has all free incidences 
in a single face. The former can be achieved as we can assume orthogonal representations to be 
tight. If the latter condition is violated by a cutvertex v, then v has two incident edges ei and ei 
and the rotation of v is in both incident faces. If both edges e\ and e2 have zero bends, we bend 
along a cycle around v in the flex graph and thus we can assume without loss of generality that e± 
has a bend. Moving v along e\ to this bend yields an orthogonal representation where v has both 
free incidences in the same face. Thus given the orthogonal representations for the blocks, we can 
simply stack them together without causing additional cost. □ 

6.3 Computing the Flow 

In the previous sections we used Tfl ow (n) as placeholder for the time necessary to compute a 
minimum-cost flow in a flow network of size n. Most minimum-cost flow algorithms do not consider 
the case of multiple sinks and sources. However, this is not a real problem as we can simply 
add a supersink connected to all sinks and a supersource connected to all sources. Unfortunately, 
the resulting flow network is no longer planar. Orlin gives a strongly polynomial time minimum- 
cost flow algorithm with running time 0{m\ogn{m + nlogn)), where n is the number of vertices 
and m the number of arcs [TH] • Since our flow network is planar (plus supersink and supersource) 
the number of arcs is linear in the number of nodes. Thus with this flow algorithm we have 
T flow (n) G 0(n 2 log 2 n). 

Cornelsen and Karrenbauer give a minimum-cost flow algorithm for planar flow networks with 
multiple sources and sinks consuming C(^/xnlog 3 n) time [3], where x is the cost of the resulting 
flow. Since the cost functions in an instance of OptimalFlexDraw may define exponentially large 
costs in the size of the input, we cannot use this flow algorithm in general to obtain a polynomial 
time algorithm. However, in practice it does not really make sense to have exponentially large costs. 
Moreover, in several interesting special cases an optimal solution has cost linear in the number of 
vertices. We obtain the following results. 

Corollary 3. A convex instance G of OptimalFlexDraw can be solved in 0(n 4 log 2 n) and 
0(i/xn 3 log 3 n) time, where \ is the cost of an optimal solution. The running time can be improved 
by a factor of 0{n) for biconnected graphs. 
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7 Conclusion 



We presented an efficient algorithm for the problem OptimalFlexDraw that can be seen as the 
optimization problem corresponding to FlexDraw. As a first step, we considered biconnected 
4-planar graphs with a fixed embedding and showed that they always admit a nice drawing, which 
implies at most three bends per edge except for a single edge on the outer face with up to four 
bends. 

Our algorithm for optimizing over all planar embeddings requires that the first bend on every 
edge does not cause any cost as the problem becomes A^P-hard otherwise. Apart from that restric- 
tion we allow the user to specify an arbitrary convex cost function independently for each edge. 
This enables the user to control the resulting drawing. For example, our algorithm can be used 
to minimize the total number of bends, neglecting the first bend of each edge. This special case 
is the natural optimization problem arising from the decision problem FlexDraw. As another 
interesting special case, one can require every edge to have at most two bends and minimize the 
number of edges having more than one bend. This enhances the algorithm by Biedl and Kant [I] 
generating drawings with at most two bends per edge with the possibility of optimization. Note 
that in both special cases the cost of an optimal solution is linear in the size of the graph, yielding 
a running time in 0(nz log 3 n) (0(n,2 log 3 n) if the graph is biconnected). 
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